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Abstract Axisymmetric resistive MHD simulations for radially self- similar ini- 
tial conditions are performed, using the NIRVANA code. The magnetic diffusivity 
could occur in outflows above an accretion disk, being transferred from the underly- 
ing disk into the disk corona by MHD turbulence (anomalous turbulent diffusivity), 
or as a result of ambipolar diffusion in partially ionized flows. We introduce, in 
addition to the classical magnetic Reynolds number Rm, which measures the im- 
portance of resistive effects in the induction equation, a new number Rb, which 
measures the importance of the resistive effects in the energy equation. We find two 
distinct regimes of solutions in our simulations. One is the low-resistivity regime, in 
which results do not differ much from ideal-MHD solutions. In the high-resistivity 
regime, results seem to show some periodicity in time-evolution, and depart signifi- 
cantly from the ideal-MHD case. Whether this departure is caused by numerical or 
physical reasons is of considerable interest for numerical simulations and theory of 
astrophysical outflows and is currently investigated. 
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1 Introduction 

In Vlahakis & Tsinganos (1998) general classes of self-consistent ideal-MHD so- 
lutions have been constructed. In Vlahakis et al. (2000) Blandford & Payne (1982) 
model was analysed, and the problem with the terminal wind solution (which was 
not causally disconnected from the disk) has been solved. The common deficiency 
of all radially self- similar models, a cut-off of the solution at small cylindrical radii 
and also at some finite height above the disk because of a strong Lorentz force close 
to the system's axis has been corrected numerically. A search in the numerical sim- 
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Fig. 1 The initial setup, which is slightly modified analytical solution, is shown in the Left panel. 
The solid lines represent logarithmically spaced isocontours of density. It is also shown in colour 
grading, in red to violet colour, for the logarithm of density -1 to -4, respectively. In the Right panel 
shown is, in the same grading, the solution with large magnetic diffusivity. It does not reach station- 
ary state, and shows some periodicity in time evolution. The dashed lines depict poloidal magnetic 
field lines, and the dotted lines depict the fast magnetosonic, Alfven and slow magnetosonic critical 
surfaces, top to bottom, respectively. 
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ulations for solutions at larger distances from the disk has been performed in Gracia 
et al. (2006) with NIRVANA code (version 2.0, Ziegler, 1998), and similar results 
were obtained also using the PLUTO code in Matsakos et al. (2008). Extension in 
the resistive-MHD has been investigated in Cemeljic et al. (2008) using the NIR- 
VANA code and some of results we present here. Results of these investigations 
could also have implications in the numerical simulations of magneto spheric inter- 
action in vicinity of the young stellar objects, where the resistivity plays important 
role. 

Our numerical simulations are initiated by the slightly modified analytical solu- 
tions for radially self- similar flow from Vlahakis et al. (2000), and then evolved in 
the resistive MHD simulations by NIRVANA code. Our initial setup is shown in the 
left panel of FigH] 



2 New characteristic number 

In addition to the magnetic Reynolds number Rm^VR/T] , which describes influence 
of the magnetic diffusivity i] in the induction equation, we introduced a new number, 
which describes the influence in the energy transport equation- see Fig|2l It can be 
written in terms of Rm and plasma beta as Rb=Rmj3/2. It is the ratio of the pressure 
term over the Joule heating term in the energy equation. When Rb is smaller or close 
to unity, which can happen even when Rm is much larger than unity, the energy 
dissipation becomes important. It might define one additional mode of resistive- 
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Fig. 2 Value of j3 /2(y7?/Vb^o) for the analytical solution used as initial setup here. This quantity 
gives the critical value of magnetic diffusivity r] that corresponds to Rb=l. 




4 Miljenko Cemeljic, Jose Gracia, Nektarios Vlahakis and Kanaris Tsinganos 

MHD solutions, indicated in our search for eventual onset of super-critical resistive 
regime. 



3 Results 



The resistive MHD jets are similar to ideal-MHD solutions for a finite range of mag- 
netic diffusivity, in which they reach a well defined stationary state. This state only 
slightly differs from the initial state, as expected, since the initial setup was slightly 
modified analytical stationary solution. Departure from the ideal-MHD regime oc- 
curs for larger values of magnetic diffusivity, above some critical value. One such 




Fig. 3 Reconnection and re- shaping of the magnetic field in the vicinity of the young stellar ob- 
ject. Initially pure stellar dipole field reshapes into the stellar and disk open field during the time- 
evolution. The time is measured in the number of rotations at the inner disk radius, which is at 
Rj-3.0 in these simulations. Without the substantial resistivity, reconnection does not occur and 
simulations terminate because of numerical reasons. 
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result is shown in the right panel of Fig[T] We note possible existence of the dis- 
tinct super-critical regime in magnetic diffusivity in our simulation for the outflows 
initialised with a self- similar analytical solutions. We also define the new character- 
istic number, which describes the influence of the resistivity on the energy transport 
equation. Physical parameters and the eventual periodicity of the super-critical resis- 
tive solutions are currently under investigation. Such solutions might be interesting 
for investigations of accretion flows in the vicinity of young stellar objects, where 
the magnetic resistivity seems to play important role. 

In the FigOwe show one case in numerical simulations of magneto spheric inter- 
action in the closest vicinity of young stellar object (Cemeljic et al., 2009). These 
simulations have been performed with code ZEUS 347, which is our resistive version 
of Zeus-3D code (Fendt & Cemeljic, 2002). Magnetic reconnection shows to play 
essential role in re-shaping the initial stellar dipole, which enables the launching of 
outflows. In our numerical simulations, if the resistivity is too small, reconnection 
will not occur. Therefore, we need to investigate parameter space for resistivity, and 
we need to understand what are the effects of large, and not only negligible or very 
small resistivity. 
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